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Abstract 

In this paper, we investigate continuous-time quantum walk on star graphs. It is 
shown that quantum central limit theorem for a continuous-time quantum walk on star 
graphs for AT-fold star power graph, which are invariant under the quantum component of 
adjacency matrix, converges to continuous-time quantum walk on K2 graphs (Complete 
graph with two vertices) and the probability of observing walk tends to the uniform 
distribution. 
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Quantum walks were introduced in the early 1990s by Aharonovich, davidovich and Zaggury 
[1]. Since then the topic has attracted considerable interest. The continuing attraction can 
be traced back to at least two reasons. First, the quantum walk is of sufficient interest in its 
own right because there are fundamental differences compared to the classical random walk. 
Next, quantum walks offer quite a number of possible applications. One of the best known is 
the link between quantum walks and quantum search algorithms which are superior to their 
classical counterparts [2, 3]. Similar to classical random walk there are two types of quantum 
walks, discrete and continuous time [1, 4]. A study of quantum walks on simple graph is well 
known in physics (for more details see [5]). Recent studies of quantum walks on more general 
graphs were described in [2, 6, 7, 8, 9, 10, 11, 12]. Some of these works studies the problem 
in the important context of algorithmic problems on graphs and suggests that quantum walks 
is a promising algorithmic technique for designing future quantum algorithms. One approach 
for investigation of continuous-time quantum walk (CTQW) on graphs is using the spectral 
distribution associated with the adjacency matrix of graphs [13, 14, 15, 16, 17, 18, 19]. Authors 
in Refs.[13, 14] have introduced a new method for calculating the probability amplitudes of 
quantum walk based on spectral distribution. In this method a canonical relation between the 
Fock space of stratification graph and set of orthogonal polynomials has been established where 
leads to obtain the probability measure (spectral distribution) adjacency matrix graph. The 
method of spectral distribution only requires simple structural data of graph and allows us to 
avoid a heavy combinational argument often necessary to obtain full description of spectrum 
of the adjacency matrix. 

In this paper we investigate CTQW on star graphs, which are discussed by Burioni et, al. 
[20, 21] as models of the Bose- Einstein condensation, by quantum probability theory. To this 
aim, first by turning the graph into metric space based on distance function, we have been able 
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to stratify and quantum decomposition [13, 23], such that the basis of Hilbert space of quantum 
walk consist of superposition of quantum kets of vertices belonging to the same stratum with 
the same coefficient (Fock space). Then by applying an isomorphism from Fock space onto 
the closed linear span of the orthogonal polynomials with respect to the measure /i we obtain 
spectral distribution adjacency matrix graph and probability amplitudes of observing walk at 
strata. Also we try to investigate quantum central limit theorem for CTQW in A^-fold star 
power graphs as N — )• oo, and show the CTQW converges to CTQW on K2 (Complete graph 
with two vertices) where the probability of observing walk tends to the uniform distribution. 
On the other hand, quantum walker induces to the graph that spectral distribution it tend to 
the Bernoulli distribution /j, — > j{Si + 5_i). 

The organization of the paper is as follows: we give a brief review of graph, star graph, 
stratification and quantum decomposition in the section 2. Section 3 is devoted to study 
CTQW for some example of star graphs via quantum probability theorey and try to investigate 
quantum central limit theorem for CTQW on their graphs. In the conclusion we summarize the 
obtained results and discuss possible development. Finally, in the appendix the determination 
of spectral distribution associated with adjacency matrix by Stieltjes transform is derived. 

2 Star graph, Stratification and Quantum decomposi- 
tion 

Let y be a non-empty set and £" be a subset of {{a, l3}\a, P E V ior a ^ The pair 
G — {V, E) is called a graph, where elements of V and E are vertices and edges of graph, 
respectively. We say that two vertices of a and /? are adjacent if {a, (3} & E and write ck ~ /3. 
A finite sequence ao, ai, is said a walk of length n if ~ a^+i for /c = 0, 1, n — 1. A 

graph is called connected if any pair of distinct vertices is connected by a walk. The degree or 
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valency of a vertex a & V is defined by 

K{a) = \{P e V\ (3 a}\ (2-1) 

where |.| denote the cardinahty. For a graph G the adjacency matrix A is given by 

1 if q; ~ 
otherwise. 

Obviously, (i) A is a symmetric (ii) elements of A take a value in {0, 1} (iii) diagonal ele- 
ments of A are 0. Conversely, for a non-empty set y, a structure graph is uniquely determined 
by a such matrix which indexed by V . On the other hand, A is considered as an operator 
acting on the Hilbert space /^(l^) in such a way that 

A\a) = Y.\(5), aey, 

where {|a)| a G forms a complete orthogonal basis of l^{y). 

In this paper, we focus on investigating CTQW on star graphs which are also discussed by 
Burioni et al. [20, 21] as models of the Bose- Einstein condensation. A star graph is obtained 
from star product of graphs which we define in the following. 

Let us consider two graphs C, v — 1,2 with adjacency matrices A" ^ each of which is 
equipped with an origin Oj, e V'^^^ . Define a matrix A as 

Aa,/?),(d,/3) = 4'i'^/3o.^4,, + ^ao/do.<j, a,«eFW, /3,/3ey(2) (2-2) 

We observe that A is symmetric, elements take a value in {0, 1} and the diagonal consists of 0. 
Therefore A is the adjacency matrix of a graph. The connected component containing (oi, 02) 
is called the star product of G*^^^ and G*^^^ and the resulted graph is a star graph [23] (see Fig.l). 
One can generalize this definition ioT v — 1,2, ...,N that we consider TV-fold star power graphs 
for our aims. The stratification is introduced [13, 23] by taking — {oi, 02, on) as the origin 
and have 

00 

V=[jVk, Vi^{ae V\ d{o, a) = k}. (2-3) 

k=0 
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Here d{a, j3) stands for the length of the shortest walk connecting a and [3. According to the 
stratification (2-3), we define a unit vector by 

1'^^) = ^ E l^>«>. (2-4) 

where |A;, a) denotes the eigenket of the a-th vertex at the stratum k and let T{G) the closed 
subspace of /^(V^) spanned by Moreover, the stratification (2-3)give rise to a quantum 

decomposition: 

A = + + (2-5) 

If r(G) is invariant under the actions of the quantum components quantum components A^, 
e e — , 0}, then there exist two Szego- Jacobi sequences {ijJk}'^=i and {cujt}^! derived from 
A, such that 

A-^\4>k)^y/^Mk+i). k>0 (2-6) 

A-|0o) = O, A-\cPk) = y/^J^\(f)k-i), k>l (2-7) 

= k>0. (2-8) 

The above coefficients (i.e., Szego- Jacobi sequences ) are obtained from geometric feature of 
graph [13, 23]. Then (T {G) , A'^ , , A°) is an interacting Fock space associated with Szego- 
Jacobi sequences {u!k,ak}. 



3 CTQW on star graph via quantum probability theory 

In the study of CTQW on graphs, the spectral distribution of A in the vacuum state |0o) plays 
an essential role [13], which is by definition a probability distribution /i uniquely specified by 

= {(t>o\A"'\<Po) = / x^fi{dx), m = 0, 1, 2, (3-9) 

where, according to [13, 14, 23], {A"^) coincides with the number of m-step walks starting and 
terminating at origin point a. For analyzing the spectral distribution // of adjacency matrix 
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A, we use the method of quantum decomposition which is a powerful tool. The spectral 
distribution is determined by applying the canonical isomorphism from the interacting Fock 
space onto the closed linear span of orthogonal polynomials determined by Szcgo- Jacobi 
sequences {ujk, ak}- In fact the determination of /i is the main problem in the spectral theory 
of operators, where in the case of star graphs is quite possible by using the Stieltjes method, as 
it is explained in appendix A. Then by using the quantum decomposition relations (2-5)-(2-8) 
and the recursion relation of polynomials P„(n)(A-l), the other matrix elements as 

{(f>k\A'^\<Po) = , ^ / x^Pk{x)i^{dx), 0,1,2,.... (3-10) 

y/UJiLU2...UJk J 

One of our main goals in this paper is the evaluation of probability amplitudes of CTQW by 
using Eq.(3-10), such that we have 

qk{t) = (</>fc|e-^*^|0o) = , ^ / e-''^Pk{x)ii{dx), (3-11) 

^yUJlUJ2...UJk J 

where |q'fe(t)P is the probability of observing the walk at the stratum k at time t. The con- 
servation of probability Z]fcl?ik(OP = 1 follows immediately from Eq. (3-11) by using the 
completeness relation of orthogonal polynomials Pn{x). In the appendix A reference [13] is 
provided the walker has the same probabihty at the all sites belonging to the same stratum, 
i.e., we have \qikit)\'^ — '^'^'\vl\ ' ^ ^ ^) where \qik{t)\'^ denotes the probability of the walker 
at the z-th vertex of fc-th stratum V^. Investigation of CTQW via spectral distribution method, 
which is introduced as a new method for calculating the probability amplitudes quantum walk 
(for more details see [13] ), allows us to avoid a heavy combinational argument often necessary 
to obtain full description of spectrum of the Hamiltonian. 

We can now investigate CTQW on star graphs. In the first, we calculate CTQW on two 
simple star graph and study quantum central limit theorem for them. In the end, we study 
this question for star lattice. 
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3.1 Examples 1. 

In this example we investigate the CTQW on the star finite path graph with three vertices as 
Fig. 2. Then we have 

ioi = N, U2 = l; ctn = (3-12) 

and 

Therefore the amphtudes probabihty of CTQW on this graph are as follows: 

qo{t) = J e-'^'fi{dx) = ^^[1 + Arcos(v/]vTTt)] 



Qiit) = 4^1 ^e-'^'Kdx) = -zy^^sin(x/iVTTt) 



q,{t) = ^J{x'- N)e-^^'^{dx) = [-1 + cos(^iVTTt)]. (3-14) 



J {N + 1)VN' 

Quantum central limit theorem 

Having studied CTQW on this graph for A*" arbitrary, we investigate CTQW in the limit 
of large N — > oo which in fact we discuss this question as a quantum central limit theorem 
for CTQW. To state a quantum central hmit theorem for CTQW we have[16] 

qk{t)= lim ((/)fc|eW|0o) = lim ^ ^ [ Pfc(x)e-^"*/^/Xjv(x)cia; (3-15) 

Then we obtain the amplitudes probability as 

f —itx r —itx X TV" 

= ^(e-** + e^*) = cost, (3-16) 

1 r -itx 1 (• -itx 1 TV 
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he-'' - e'') ^ -isint, (3-17) 

q2{t)=0, (3-18) 

where the probabihty of the observing walk (|gj(t)p, for i — 0, 1 ) reaches the uniform distri- 
bution periodically at t = kn + j for k & Z . 

Indeed, in the limit of large N — )• oo(e.i., quantum central hmit theorem ), the CTQW on 
this graph is reduced to CTQW on K2 (Complete graph with two vertices). In other words, 
the spectral distribution fx^ associated with the adjacency matrix of the graph tend to the 
Bernoulli distribution (/xjv — > |((^i + (^-1)), where is the spectral distribution of K2 graphs. 

3.2 Examples 2. 

In the second example we consider CTQW on the star graph as Fig.3. Then we have 

LOi = 2N, cua = 2; = (3-19) 

and 

a;2 _ 2 1 N f \ 

The amplitudes probability of CTQW on this graph are as: 

q,{t) = j e-V(c?x) = ^^[1 + iVcos(v^2(iV+l)t)] 



q,{t) = ^J xe-'^'tx{dx) = -^^^sin(y2(ivTT)^) 

q2it) = I (x'^ - N)e-'''' nidx) = — — =[-1 + 2cos(^/NTTt)]. (3-21) 

Also in this example, quantum central limit theorem for CTQW is reduced on CTQW on 
complete graph K2. 
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3.3 Examples 3. 
Star Lattice 

The star lattice is the TV-fold star power of finite path graph, see Fig.4. One can show that 
the two sequence {u>k} and {cufe} obtain as 

cui^N, a;2 = a;3 = = 0:^ = (3-22) 

Substituting sequence Uk and ak in (5.22), the Stieltjes transform Gnj^{z) of spectral distribu- 
tion jjLjsi takes the following form 

G,A^) = \ (3-23) 

z 1 

1 - 



In order to evaluate the continued fraction, we need first to evaluate the following infinite 
continued fraction defined as 

G,{z) = , ^ 1 = i— , (3-24) 

where by solving above equation we have 

G,{z) = ^{z-^/^^). (3-25) 

Then by substituting (3-25) in (3-31), we obtain the following expression for the Stieltjes 
transform of fj, 

G (z) ^ _ H2-N)z-NV^ 

Finally by applying Stieltjes inversion formula, we acquire the absolutely continuous part of 
spectral distribution /xjv as follows 



1 TV a/4 - 1-2 
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Now we investigate CTQW on some of the known infinite graphs can be obtained from star 
lattice by appropriate choice of A^" as: 

A. For N — 1 we obtain finite path graph. Then we have 

I2^(x) = - ; -2<x<2, (3-28) 

where the amphtudes probabihty for CTQW on this graph are 

Qkit) - t'iMt) + Jk+2{t)) = 2i\k + 1):^^; - 0, 1, ... (3-29) 

(for more detail see Ref.[13]). 

B. For N — 2 we obtain infinite line graph where the two sequences, spectral distribution 
and the amplitude probability for observing walk at strata at time t are 

a;i = 2, = a;3 = • • • = 1; afe = 

ii^ix) = -^=i= ; -2 < X < 2, 

qo{t)^Mty, qk{t)^i''V2Jk{t), k>l, (3-30) 

(for more detail see Ref. [13]). 

C. For A'" > 3 the total mass of Hn{x) is less than one (i.e., / iJiN{x)dx < 1). In fact, it 
contains a discrete measure where by a simple computation one can obtain as [24] 



1 N\/A - x^ N -2 

""(^^ = 2i N^-(N-l).^ + 2F3^(** + • -2 ^ ^ ^ 

Quantum central limit theorem 

Now we consider behavior CTQW on this graph for large (i.e., N — >■ oo) where we 
discuss this question as a quantum central limit theorem for CTQW. 
The quantum central limit theorem for CTQW is 

Qkit) = hm {cf>,\e-'^'/^\4>o) = lim ^ / P,{x)e-'^'/^ ii^ix)dx 

N — >oo N — >oo y JS J 
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= hm ^= / Pk(x)e — — — — + — (5 N +S -N )]dx 

= lim ^ / Pfc(a;)e-^^*/^ (5 n + 5 -n )dx 

iv-.oo 2N -2\ ^ ViV j ^ ' 

By a simple computation one can show that hmjv >oo — = for A; > 2. Then we have 

9oW = ^(e-'* + e^*)=cosi, 

Qkit) =0 for k>2. (3-33) 

We see that quantum central hmit theorem for CTQW on the star lattice is reduced on complete 
graph K2 and the amplitudes probability reach the uniform distribution. Really, for every star 
graph that is invariant under the quantum component A^, e & {+,—,0}, one can show that 
quantum central hmit theorem for CTQW on it is reduced CTQW on complete graph K2. 

4 Conclusion 

In this work, we have investigated CTQW on star graphs, which are discussed by Burioni et, 
al. [20, 21] as models of the Rose- Einstein condensation, by quantum probability theory. First, 
we have equipped graph to metric based on distance function. Then we have stratified graph 
and according to stratification accomplished quantum decomposition of adjacency matrix. The 
basis of Hilbert space of quantum walk form a Fock space where by applying an isomorphism 
from it onto the closed linear span of the orthogonal polynomials we have obtained spectral 
distribution adjacency matrix graph and probability amplitudes of observing walk at strata. 
Also, we have studied quantum central limit theorem for CTQW on star graphs and shown the 
CTQW converges to CTQW on K2 (complete graph with two vertices), such that the prob- 
ability of observing walk tend to the uniform distribution. Indeed, we can generalize CTQW 
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on star graphs which are invariant under the quantum component A^, e e {+, — , o}, quantum 
central hmit theorem for them induce in CTQW on graphs that the spectral distribution it 
tend to the Bernoulli distribution. 



Appendix A 

Determination of spectral distribution by the Stieltjes transform 

In this appendix we explain how we can determine spectral distribution fj,{x) of the graphs, 
by using the Szego-Jacobi sequences {{uk},{ak}). To this aim we may apply the canonical 
isomorphism from the interacting Fock space onto the closed linear span of the orthogonal 
polynomials determined by the Szego-Jacobi sequences ({a;^}, {oj}). More precisely, the spec- 
tral distribution n under question is characterized by the property of orthogonalizing the 
polynomials {Pn} defined recurrently by 

Po{x)^l, Pi{x)^x-ai, 

xPn{x) = Pn+i{x) + an+iPn{x) + a;„P„_i(x), (A-1) 

for n > 1. 

As it is shown in [25] , the spectral distribution can be determined by the following identity: 

G (z)^ [ = - = ^"^-'^^^ = y (A-2) 

" Jr z-x z-ai- , -^-2 Pniz) fr[z-xi' 



where Gfj,{z) is called the Stieltjes transform and Ai is the coefficient in the Gauss quadrature 
formula corresponding to the roots xi of polynomial Pn{x) and where the polynomials {Qn ■*} 
are defined recurrently as 

Qi\x) = X - a2, 

xQ(^\x) = Qi^Ux) + a,,+2Q^^^x)+u;^+^Qi^Ux), 
for n> 1. 
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Now if G^{z) is known, then the spectral distribution jji can be recovered from Gn{z) by 
means of the Stieltjes inversion formula: 



Substituting the right hand side of (A-2) in (A-3), the spectral distribution can be determined 
in terms of xi,l — 1,2,..., the roots of the polynomial Pn{x), and Guass quadrature constant 
A;, 1 — 1,2, ... as 



( for more details see Refs. [13, 14, 25, 26].) 
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Figure Captions 
Figure-1: Star graph. 
Figure-2: Figure. 2 
Figure-3: Figure. 3. 
Figure-4: Star lattice for N 



